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We report the first direct observation of the oxygen-isotope ( le O/ 18 0) effect on the in-plane 
penetration depth A a ;, in a nearly optimally doped YBa2Cu307_i film using the novel low-energy 
muon-spin rotation technique. Spin polarized low energy muons are implanted in the film at a 
known depth z beneath the surface and precess in the local magnetic field B(z). This feature allows 
us to measure directly the profile B(z) of the magnetic field inside the superconducting film in the 
Meissner state and to make a model independent determination of \ a b- A substantial isotope shift 
A\ a b/Xab = 2.8(7)% at 4 K is observed, implying that the in-plane effective supercarrier mass m* ab 
is oxygen-isotope dependent with Am*i/m^ = 5.5(1.4)%. 

PACS numbers: 76.75.-K, 74.72.Bk, 74.78. Bz, 82.20.Tr 



One of the fundamental questions concerning the 
physics of cuprate high-temperature superconductors 
(HTSC) is whether the electron-phonon interaction plays 
an essential role in these systems or not. The high val- 
ues of the superconducting transition temperature T c and 
the observation of only a tiny oxygen-isotope shift of 
T c in optimally doped HTSC (see, e.g. Q) were taken 
as important arguments to propose alternative pairing 
mechanisms of purely electronic origin. The conventional 
phonon-mediated theory of superconductivity is based on 
the Migdal adiabatic approximation in which the effective 
supercarrier mass m* is independent of the mass M of 
the lattice atoms. However, if the interaction between the 
carriers and the lattice is strong enough, the Migdal adia- 
batic approximation breaks down and m* depends on M 
(see, e.g. Q). The ideal experiment to explore a possible 
coupling of the supercarriers to the lattice is an isotope 
effect study of the magnetic field penetration depth. For 
HTSC, which are superconductors in the clean limit, the 
in-plane penetration depth A a f, obeys the relation: 



!Aab <* n 8 /m* ab , 



(1) 



where n s is the superconducting charge carrier density 



and 



is the in-plane effective mass of the charge 



carriers. Note that there is no such simple relation 
for T c . Since n s was found to be predominantly iso- 
tope independent H S S H 11 @, isotope ef- 
fect experiments on the penetration depth turned out 
to be a unique tool to investigate unconventional lat- 
tice effects in HTSC. Previous oxygen-isotope effect 
studies of the penetration depth in YBa2Cu307_,5 
La 2 -zSr a Cu0 4 0,0,11, Bi 1 .«PbQ. 4 Sr 2 Ca 2 Cu 3 Oio+,5 
and Yi_ 2; Pr. c Ba2Cu307_ ( 5 M showed indeed a pro- 
nounced oxygen isotope dependence of the supercarrier 
mass. In all these experiments, however, the penetra- 
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tion depth was determined indirectly from magnetization 
measurements [3J , the Meissner fraction [4J, |5j , magnetic 
torque [6J, and muon-spin rotation (/iSR) experiments 
0. The expression "indirectly" means that all these 
techniques require a model to extract the absolute value 
of the penetration depth A. In contrast, measurements 
of the exponential field decay (B(z) = Bo exp(— z/X)) in 
the Meissner state allow to determine A in a direct way. 
Recently, such a direct measurement of the in-plane mag- 
netic penetration depth A a f, in a thin YBa2Cu307_5 film 
was performed 10] at the Paul Scherrer Institute (PSI, 
Switzerland) by using the novel low-energy /iSR (LE/xSR) 
technique . In this experiment the muon implantation 
depth was controlled by the variation of the incoming 
muon energy and the magnitude of the field was mon- 
itored by the muon spin precession frequency, alike in 
standard /iSR. 

In this Letter, we report on the oxygen-isotope 
( 16 0/ 18 0) effect (OIE) on the in-plane magnetic field 
penetration depth in a high-quality YBa2Cu307_,5 film 
near optimal doping measured directly by LE/xSR. The 
experiments revealed that in the Meissner state the mag- 
netic field B(z) in the 16 substituted sample decreases 
stronger with distance z than in the 18 O substituted sam- 
ple, clearly demonstrating that 16 A a b < 18 X a b- 

The samples used for the LE/iSR experiments were 
high quality YBa2Cu3 07_a films with an area of 2x3 cm 2 
and a thickness of 600 nm supplied by Theva 12]. The 
films were grown by thermal coevaporation of the con- 
stituents onto a single-crystal barium titanate substrate. 
The oxygen-isotope exchange was performed by anneal- 
ing the samples in 18 02 gas. In order to ensure that the 
16 and ls O substituted samples are subject of the same 
thermal history, the annealing (in 16 02 and 18 02) was 
performed simultaneously. The exchange process was 
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carried out at 600°C during 25 h, followed by slow cool- 
ing (20°C/h) in order to oxidize them completely. The 
fraction of the 18 content in the film was estimated to 
be 95%, corresponding to the 18 content in the initial 
gas used for the annealing process. The oxygen content 
(7-8) and the quality of the films were estimated from 
lmT field-cooled (FC) SQUID magnetization measure- 
ments. The observed T c onsets are 16 T C = 89.3(1) K 



and 18 T C = 89.1(1) K. The transition widths are ~0.9 K 
for both samples. A comparison of the T c 's with the 
known experimental T c vs 8 curves for 16 and 18 sub- 
stituted YBa 2 Cu 3 7 _ 5 .13] yields 8 = 0.150(5) for both 
films. While 8 indicates that the samples are slightly un- 
derdoped the narrow transition shows that the oxygen 
distribution over the films is quite homogeneous. The 
results of the OIE on T c are summarized in Table [I] 



TABLE I: Summary of the OIE results for YBa2Cu307_i extracted from the experimental data (see text for an explanation). 



Method 


Sample 


IBj. IHj, 

(K/ (K)" 


AT C /T C A\ ab /\ ab 
(%) (%) 


Am* b /m* ab 
(%) 


LE/iSR 


Thin film 


89.3(1) 89.1(1) 


-0.22(16) 2.8(7) 


5.5(1.4) 


Magnetization 


Fine powder 


91.66(3) 91.42(3) 
91.71(3)" 91.45(1)' 


-0.26(5) 3.0(1.1) 
-0.28(5) 2.4(1.0) 


6.0(2.2) 
4.8(2.0) 



"results for the back-exchange ( ls O^ ls O) sample 
''results for the back-exchange ( lci O^ ls O) sample 



The transverse-field LE/iSR experiments were per- 
formed on the 7rE3 muon beam line at PSI. A weak 
external magnetic field Bq = 9.2 mT was applied par- 
allel to the sample surface after the sample was cooled 
in zero magnetic field from a temperature above T c to 
4 K. Spin-polarized muons were implanted at a depth 
ranging from 20-150 nm beneath the surface of the film 
by varying the energy of the incident muons from 3 to 
30 keV. For each implantation energy a time-differential 
/iSR spectrum was measured. The muon implantation 
depth profile n(z) for the given implantation e nerg y was 
calculated using a Monte-Carlo code TRIM.SP [3|. The 
reliability of the calculated n(z) was crosschecked by pre- 
vious LE/iSR experiments on thin metal layers [l5j. 

The experimental data were analyzed in the following 
way. For each implantation energy the average value of 
the magnetic field B and correspondent average value 
of the stopping distance z were extracted. The value of 
B was taken from the fit of the time evolution of the 
muon-spin polarization spectrum by using the Gaussian 
relaxation function: 



a(t) = a exp(-cr 2 £ 2 /2) cos(^Bt + . 



(2) 



where ag is the initial asymmetry, a is the Gaussian relax- 
ation rate, and 7 = 27r • 135MHz/T is the gyromagnetic 
ratio of the muon. Note that at our level of statistics 
(~ 5 x 10 5 events per spectrum) the fit of experimental 
data with Eq. J2J) satisfies the \ 2 criterium. The value 
of z was taken as the first moment of the emulated n(z) 
distribution. Results of this analysis for the 16 O and 



18 substituted YBa2Cu307_5 films are shown in Fig.^ 
The data points for the 18 O film are systematically higher 
than those for the 16 O film, showing that 16 A a fc < 18 A a fc. 
The solid lines represent a fit to the B data by the func- 
tion: 



B(z) = B a 



cosh[(t - z)/X ab ] 
cosh(t/X ab ) 



(3) 



This is the form of the classical exponential field de- 
cay in the Meissner state B(z) = B cxp(—z/X ab ), mod- 
ified for a film with thickness 2t with flux penetrat- 
ing from both sides. The value of z was corrected by 
zq = 8 nm, corresponding to a "dead layer", which 
mainly arises from the surface roughness of the film |lp| . 
Fits with Eq. © to the extracted w B(z) and ls B(z) 
yield 16 A o6 (4K) = 151.85(75) nm and 18 A ab (4K) = 
155.82(68) nm. Taking into account a 18 content of 
95%, the relative shift was found to be AX ab /X ab = 



( 18 X ab 



16 



Xab)/ w X ab = 2.8(7)% at 4 K. This value 



is consistent with previous results for optimally doped 
YBa 2 Cu 3 7 -5 @ and Bii. 6 Pbo.4Sr 2 Ca2Cu3Oio+<s0 ob- 
tained from indirect magnetization measurements. Here 
we would like to mention two important points: (i) The 
good agreement between the experimental points and 
theoretical curves (Eq. |J3J) indicates that the films are 
homogenous over the whole thickness, (ii) X ab obtained 
from LE^SR is model independent. It is extracted from 
the measured field profile in the Meissner state using only 
one fit parameter A a h. 

We also performed additional OIE experiments on 
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FIG. 1: Magnetic field penetration profiles B(z) for a le O sub- 
stituted (closed symbols) and a 18 O substituted (open sym- 
bols) YBa2Cu307_<5 film measured in the Meissner state at 
4 K and an external field of 9.2 mT, applied parallel to the 
surface of the film. The data are shown for implantation en- 
ergies 3, 6, 10, 16, 22, and 29 keV starting from the surface 
of the sample. Solid curves are best fits by Eq.JHJ. 



X a b based on measurements of the Meissner fraction in 
fine YBa2Cu307_,5 powder. The powder samples were 
ground for about 60 min and then passed through a 
10 /jm sieve. The oxygen exchange procedure was anal- 
ogous to the one described above for the thin film sam- 
ples. The 18 content in the sample, as determined 
from the change of the sample weight after the iso- 
tope exchange, was found to be 89(2)%. The value of 
the Meissner fraction / was calculated from 1 mT FC 
SQUID magnetization measurements. The absence of 
weak links between grains was confirmed by the linear 
magnetic field dependence of the FC magnetization mea- 
sured at 5 K in 0.5 mT, 1 mT, and 1.5 mT. The tem- 
perature dependence of the magnetic field penetration 
depth A G ff (powder average) was calculated from / un- 
der the assumption that the sample grains are spherical 
[lif . The in-plane penetration depth X a b(T) (Fig. was 
determined from the measured A c ff (T) using the relation 
A ff = 1.31A a f, which holds for highly anisotropic super- 
conductors (Xc/Xab > 5) 17]. Taking into account a 18 
content of 89%, the relative shift at 4 K was found to 
be AA a h/A a h = 3.0(1.1)% in a good agreement with the 
LE/iSR data (see Table HJ. I n order to substantiate the 
intrinsic character of the observed OIE on X a b, we per- 
formed a back exchange experiment. For this purpose the 
16 sample was annealed in 18 02 gas ( 16 0^ 18 0) and the 
18 sample in 16 C>2 gas ( 18 0^ 16 0). The annealing pro- 
cedure was analogous to the one described above for the 



thin film samples. The results of the back exchange ex- 
periments are also shown in Fig. [2] (cross symbols). All 
the results of the OIE on T c and X a b are summarized in 
Table H| 

According to Eq. (JTJ the OIE on X a b is due to a shift 
in n s and/or m* ab : 



AA afe 

Xab 



1 Am 
2~m 



a h 



1 An s 

2 n. 



(4) 



Previous OIE studies in HTSC's 0, 0, U 0, 0, @ 
showed a pronounced OIE on m* ab with negligible OIE 
on n s . In the present work we provide further evi- 
dence of this scenario from measurements of the nuclear 
quadrupole resonance (NQR) frequency (vq) of the plane 
(v Q l ) and the chain (v Q h ) 63 Cu in 16 O and 18 O substituted 
YBa2Cu307_,5 powder samples. It is known |18lll9] that 
in HTSC Vq and Vq K are very sensitive to changes in 
the density of mobile carriers n in the superconduct- 
ing Cu02 planes (n pl ) and CuO chains (n ch ). Both Uq 
and Vq 1 increase with doped hole concentration roughly 
20 MHz per one doped hole per Cu atom 20]. The NQR 
measurements of 63 Cu at 94 K show that within 5 kHz 
error bar limits the plane copper and the chain copper 
NQR frequencies are the same in 16 O and 18 O substituted 
YBa 2 Cu307_«5 powder samples (see Table HTf . This im- 
plies that by the oxygen substitution n pl and n ch change 
for less than 3 x 10 -4 hole per Cu atom. In optimally 
doped YBa2Cu307_,5 one has two plane and one chain 
Cu atom in the unit cell. From the NQR measurements 
we can thus conclude that the change of the hole number 
per unit cell is less than 3 ■ 3 x 10~ 4 ~ 10~ 3 hole. Since 
in optimally doped YBa2Cu307_,s there is approximately 
one doped hole per unit cell, the relative change of hole 
density An/n at the oxygen exchange must be less than 
10~ 3 . The NQR results are summarized in Table ITT1 For 
YBa2Cu307_5 it was shown [il| that n s « n at T < 
70 K. Based on this observation we thus can conclude 
that |l/2 • An s /n s | < 0.05% in Eq.Q. Consequently, 
the main contribution to the OIE on X a b has to come 
from the isotope dependence of the in-plane charge car- 
rier mass m* h , so that AX a b/X a b ~ (Am* b /m* b )/2. With 
AA afc /A afc = 2.8(7)% we obtain Am* ab /m* ab = 5.5(1.4)% 
at 4 K (see Table HJ. This result is remarkable in spite of 
the fact that the observed OIE on T c in optimally doped 
YBa 2 Cu307_«5 is rather small (AT C /T C = -0.26(5)%, see 
Tabled). 

Note that such an isotope effect on m* ab is not expected 
for a conventional weak-coupling phonon-mediated BCS 
superconductor. In fact in HTSC the charge carriers are 
strongly coupled to optical phonons, as demonstrated by 
measurements of the static and high-frequency dielectric 
constants [13] , photoemission and tunnelling [24| ex- 
periments. Strong interaction between the charge carri- 
ers and the lattice ions leads to a break down of the adi- 
abatic Migdal approximation 0], and consequently the 
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TABLE II: Results of plane and chain 63 Cu NQR in 16 and 
18 substituted YBa2Cu307_a powder samples at 94 K (see 
text for an explanation). 



lfe I8 Q 

vq Linewidth vq Linewidth Avq/vq 

(MHz) (KHz) (MHz) (KHz) (%) 

Planes 31.580(5) 336(10) 31.580(5) 340(10) 0.00(2) 

Chains 22.050(5) 161(5) 22.050(5) 161(5) 0.00(3) 



coupled to the lattice. 
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effective supercarrier mass m* depends on the mass of 
the lattice atoms. To our knowledge there are just a few 
theoretical models which predict an OIE on the effective 
carrier mass m* (see e.g. Refs. 13. Eil Erj E7I) . 
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2: Temperature dependence of X a 
(0) for le O and ls O substituted YBa2Cu 3 7 -,5 fine pow- 
der samples as obtained from low-field SQUID magnetization 
measurements. The inset shows the low-temperature region 
between 10 K and 60 K. The reproducibility of the oxygen ex- 
change procedure was checked by the backexchange (crosses). 



In conclusion, we used LE/iSR to measure directly the 
oxygen isotope ( 16 0/ 18 0) effect on the in-plane mag- 
netic field penetration depth A a b in optimally doped 
YBa2Cu307_,5 films. The OIE on A a & at 4 K was found 
to be AX a b/Xab = 2.8(7)%. The intrinsic character of 
the OIE on X a b was confirmed by back exchange ex- 
periments on fine powders with low-field SQUID mag- 
netization measurements. It is concluded that the OIE 
arises mainly from the isotope dependence of the in-plane 
charge carrier mass m* ab with Am* b /m* fc = 5.5(1.4)% at 
4 K. This finding implies that even in optimally doped 
cuprate superconductors for which only a small isotope 
effect on T c is observed, the supercarriers are strongly 



